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Abstract. The Monte Carlo event generator PHOKHARA, which simulates hadron and muon production 
at electron-positron colliders through radiative return, has been extended to final states with three pions. 
A model for the form factor based on generalized vector dominance has been employed, which is consistent 
with presently available experimental observations. 


1 Introduction 

Measurements of form factors and cross sections for elec¬ 
tron positron annihilation in the low energy region provide 
important information on hadron dynamics. At the same 
time, they are necessary ingredients in dispersion rela¬ 
tions, which are used to predict hadronic contributions to 
the momentum dependent electromagnetic coupling and 
the anomalous magnetic moment of the muon. The mea¬ 
surements are traditionally performed by tuning the cen¬ 
ter of mass energy of an electron positron collider to the 
point of interest. As an alternative, it has been advocated 
to use the method of the radiative return m at high 
luminosity (j)- and H-meson factories. In this second case 
the collider energy remains fixed, while the invariant 
mass of the hadronic system, can be varied by considering 
events, where one or several photons have been radiated. 
For a detailed and precise analysis initial and final state 
radiation must be included and a proper description of 
the various exclusive final states through appropriate form 
factors is required. All these ingredients are contained in 


the most recent version of the Monte Garlo event gener¬ 
ator PHOKHARA4.0 [HlllI, which is based in particular 
on the virtual corrections described in iniini and which at 
present simulates production of tt+tt”, four pions 

(27r+27r“ and 7r+7r“27r^), pp, and nn |j| |I| |^||d||7| |8||^ . In 
the present work the production of three pions is consid¬ 
ered on the basis of form factors, which include already 
available information on the production cross section and 
on differential distributions in two-pion subsystems. The 
model implements three-pion production through cj, 0 and 
their radial excitations and the subsequent decay of these 
resonances into pir, p'lr and p"7v. A small isospin-violating 
component 7 * ^ uo{-^ 7 r+ 7 r“) 7 r^ is added, which is needed 
to properly describe the data. The total cross section and 
the distributions are well reproduced within this model. 
It is furthermore demonstrated that the couplings intro¬ 
duced and adopted for this purpose lead to a satisfactory 
description of T( 7 r^ ^ 77 ), of the slope parameter of the 
77 * amplitude and of the radiative vector meson de¬ 
but is in conflict with uj 


TT 


cays p 


7r^7 


7r%, 


7r^7. 
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2 A phenomenological description of 
three-pion production 

The amplitude for three-pion production through the elec¬ 
tromagnetic current is restricted by current conservation 
and negative parity to the form 

= ei^ap^q+q-qo F3^{q+,q-,qo) . ( 1 ) 
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G-parity dictates dominance of the isospin-zero com¬ 
ponent of the electromagnetic current , which will be dis¬ 
cussed in a first step. The small isospin-one admixture will 
be discussed subsequently. The form factor is con¬ 

structed under the assumption that the virtual photon 
couples to the uo- and 0-meson, whose subsequent transi¬ 
tion to three pions is dominated by the p{^ 2 'k')'k chain 
(FigGj pnj . Taking into account radial excitations of a;, 
0, p one arrives at the form factor 


FL° (9+,9-,9o) = 

^ aij ■ BWv,{Q^) ■ (Ql, Ql,Ql) , ( 2 ) 

hj 

where Vi stands for either uj- or 0-resonances, and pj rep¬ 
resents contributions from p-mesons. From the PDG on it 



Fig. 1. Diagrams contributing to the 3-pion current: / = 0 
component. 


is clear that all uj- and 0-resonances couple to the ground 
state of the p-meson, whereas there is no indication about 
couplings to the higher radial excitations (p', p", • * * )• This 
missing piece of information can, however, be obtained to 
large extent from the known e+e“ ^ tt+tt^tt^ cross sec¬ 
tion, as shown below. 

For the function Hp we shall adopt the ansatz 


where Qf = {qj + and rrij = m^j , with i,j, /c = 

0, ±. We use propagators with constant widths for cj’s 
and 0, and energy dependent widths for p- resonances as 
predicted by P-wave p ^ tttt decays: 


m 


2 r 


^piQi^ 5 ) — Fp ^2 


Qi - 


2 13/2 




(6) 


The couplings aij are taken as real constants and we as¬ 
sume that the isospin symmetry is violated in this com¬ 
ponent only by the tt^ — tt^ mass difference. 



Fig. 2. Diagram contributing to the / = 1 component of the 
three-pion current. 


The small isospin violating amplitude is mediated by 
the 1=1 component of the electromagnetic current and 
is based on the 47r current of Refs. pmril i.e. we take the 
p — 7 and piruj couplings from the 47r current and replace 
the UJ ^ Stt transition by the isospin violating uj ^ 27t 
decay as shown in FigEl This leads to the following ansatz 

Fl=\q+,q., qo) = ■ BWUQ^)lml 

[BWp{Ql)/ml + aSVy. (Qg)/m2„] , (7) 

where 

1 55 

G^ = -^ 12.924 GeV^ 0.266 (8) 

V 2 


Hp{Q\M''-Ml)=BWp{Ql)FBWp{Q\)FBWp{Q'^_) , 

( 3 ) 

with 


Qo — Q± — {Qt + (4) 


and the Breit-Wigner form factors are 


BWviQ^) 

BWpiQ^) 


91 

ml 

Qi 


mi 


-1 + i 


.Fv 


my 


-iFi 


■ v^iFpiQj ,mj,mk) 


m„ 


,( 5 ) 


and fhp = 0.77609 GeV, Bp = 0.14446 GeV, fhpn = 1.7 
GeV, tpu = 0.26 GeV, a = —0.1, where the parameters 
are taken directly from m- At the present level of ex¬ 
perimental accuracy it is not clear whether the p" term is 
necessary for the description of the Stt current (see below). 
However, as it is a prediction coming from the An current 
we consider its contribution also here. 

The coupling can be extracted from the decay 

rate F{uj ^ tttt) (see Eq.(^2|). Using the world average 
value from the PDG HD one gets = 0.185(15). The 
total form factor is of course given by the coherent sum 

FsT,{q+, q-,qo) = Fi=°{q+, q-, qo) + Fiy{q+, q-,qo) ■ (9) 

Data on cr(e+e“ ^ tt+tt^tt^), from energy scan experi¬ 
ments [CTll5lllf)lll7lll8j consist of 217 data points, covering 
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Fig. 3. e+e TT^TT 7r° cross section obtained with fitted parameters (solid line, see text for details) vs. experimental data. 



Vs(M^y) y/s{MeV) 


Fig. 4. e^e 


TT^TT 7r° cross section obtained with fitted parameters (solid line, see text for details) vs. experimental data. 


the energy range from 660 MeV to 2400 MeV. Moreover, 
by using the radiative return method, BaBar JHI has ob¬ 
tained additional 78 data points, covering the energy from 
1.06 GeV up to almost 3 GeV. Data from the DM2 col¬ 
laboration m being inconsistent with the more accurate 
BaBar data, were not used in the fit. 


To fit the experimental data we include the following 
resonances: uj{782), uj' = a;(1420), uj" = a;(1650), 0(1020), 
p(770), p' = p(1450) and p" = p{1700). The strategy was 
to minimize the number of contributions in Eq. 0 and 
arrive at a good fit in the same time. Our best fit for the 


isospin zero component reads 
-p3r°(9+,9-,9o) = -ffp(770)(QbQ-’*5o) 

• A ■ BW^(782) (Q") + B ■ BW^(^o20) (Q") 

+ C ■ BW^(u20) (Q") + B ■ BIV.(1650) (Q") 

+ E • -Bvy,(io2o) (Q^) ■ -ffp(i45o) {Q +) Q-) Qo) 

+ F ■ BW(^(1650)(Q^) • i?p(1700)(Q+,Q-,Qo) ; (10) 

with the couplings and masses given in Tabled The errors 
are the parabolic errors calculated by the MINUIT pro¬ 
cessor MINOS and correspond to the change of Ax^ = 1. 
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Vs(Mey) 

Fig. 5. e+e“ —^ TT^TT 7r° cross section obtained with fitted 
parameters (solid line, see text for details) vs. experimental 
data. 


The value x^/d.o.f = 1.14 is at the edge of the 95% con¬ 
fidence interval. It is to large extent a result of the fact 
that the data of different experiments are only marginally 
mutually consistent as evident from Figures IHHSl Lack of 
published numerical information on differential distribu¬ 
tions does not allow for further refinements of the model. 
In particular, indications for additional contributions from 
higher radial excitations in the region of large cannot 
be substantiated by more detailed fits. Access to distribu¬ 
tions in invariant masses of pion pairs and to pion angular 
distributions would allow to study these effects. For the 
moment we set the masses and widths of the p mesons to 
the values collected in Table [3 and assume equal masses 
and widths for the neutral and charged p mesons. Never¬ 
theless qualitative comparisons of two-pion invariant mass 
distributions can be performed by superimposing the ex¬ 
perimental and the model (data from Fig. 15 of JHI and 
differential cross sections generated with PHOKHARA 
5.0) distributions (see Fig IS)). The / = 0 component alone 
predicts identical distributions in and M^±^o. The 


Table 1. Values of the couplings masses and widths obtained 
in the fit; couplings A — F \n GeV“^ masses and widths in 
MeV (see text for details). 


—I-1-1-1-1-1-1- ! i__L. ! ,.! TUT 
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MMGeV) 


Fig. 6. Two pion invariant mass distributions for four differ¬ 
ent ranges of invariant mass. The BaBar data points, 

given as events/bin, are superimposed on plots obtained by 
PHOKHARA (see text for details). 


^a;(782) 

782.4(4) 

A 

18.20(8) 

Ta;(782) 

8.69(7) 

B 

-0.87(5) small 1=1 component is concentrated in the spike at 

^0(1020) 

1019.24(3) 

C 

-0.77(5) = m^. This channel starts to contribute for 

F</)(1020) 

4.14(5) 

D 

-1.12(4) values above 1 GeV only. The model seems to describe 

^c.u(1420) 

1375(1) 

E 

-0.72(10) the distributions reasonably well. However the following 

TL;(1420) 

250(5) 

F 

-0.59(4) deviations are observed: In the lowest region (0.75 GeV 

^c.u(1650) 

1631(6) 


< <0.82 GeV) the experimental results for the dis- 

Tkcieso) 

245(13) 

X^/d.o.f 

1.14 tributions in M^+^- and respectively, seem to dif- 


fer in the upper range, an indication of isospin-violation, 
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Table 2. The masses and widths of p resonances used in the 



p(770) 

p(1450) 

p(1700) 

m (GeV) 

0.77609 

1.465 

1.7 

r (GeV) 

0.14446 

0.31 

0.235 


that cannot be reproduced by our ansatz. In the large 
range (1.4 GeV < < 1.8 GeV) an excess is observed 

in both charge modes for masses of the two-pion system 
between 1 GeV and 1.2 GeV. A similar excess is not ob¬ 
served in the pion form factor. 



Fig. 7. Diagrams contributing to 7 r° ^ 77 decay. 


3 Meson couplings and partial decay widths. 


From the results of the fit we can evaluate the meson cou¬ 
plings separately, combining the following relations 


A — • pujirp ’ Qp'K'K 

B ’ 9(j)7Tp ' 9p7T7T 

C — • g^j'T^p ’ gp'K'K 

D = 2g^"j • guj'fq^p • gp^n 
E — ^9(I>7 ’ 9(f>7Tp' ’ 9p'7T7T 

F = 2g^//j • guj"'Kp" * gp"7r7r 5 


with information about partial decay widths. From the 
known decay widths of p^^uj and <p to two pions one de¬ 
termines the couplings gvirir {V = p^^uo^cj)): 


Fv 


_ 2 


1 - 


4m^ , 


-i 3/2 


m 


V J 


( 12 ) 


Similarly, the gy^y couplings can be found from the mea¬ 
sured values of F(V e+e“): 




V^e+e - 9v 


^ Srriy 


(13) 


For the u' and u" we do not know the partial decay widths 
and the couplings cannot be extracted separately. The nu¬ 
merical values of the couplings obtained from Eq. and 
partial decay widths [U are collected in Table 0 


Table 3. Values of the three- and two-particle couplings; 
gvmr is dimensionless, gy^ in GeV^, gvirp in GeV“^. 


PpTTTT 

5.997(32) 

9pi 

0.1212(13) 

PcUTTTT 

0.185(15) 

guj-y 

0.03591(37) 

P^TTTT 

0.0072(6) 

94)1 

0.0777(7) 

goonp 

42.3(5) 

9u}'i ‘ 9u)'7rp 

-0.064(8) 

gcpTvp 

-0.93(5) 

9uj"i ■ 9uj"7rp 

-0.093(3) 


Having extracted the couplings, we are able to predict 
many physical quantities that have been measured already 
and check the model. In particular we will investigate if 


various meson-photon interactions can be modeled, as in 
the Stt case, by three-meson couplings and vector meson- 
photon mixings with the couplings and propagators as in¬ 
troduced above (for a review of alternative models see pui 
EH). As one can recognize our model is an extension to 
higher radial excitations of the model outlined in m- 
Let us start with the decay width of tt^ ^ 77 . The 
model is based on the diagrams shown in FiglU As in¬ 
dicated by the fit, contributions from p' — 7 and p" — 7 
mixings are not required at the present level of precision. 


Table 4. Mean life time for tt® ^ 27 in seconds and decay 
rates for ^ 7 r °7 in MeV as obtained within our model 

compared to experimental results uu. 



model 

experiment 

r{7v'^ 77) 

6.6(3)-10“"' 

8.3(6)-10""'' 

r(p" ^ Try) 

0.078(3) 

0.090(20) 

r{(jJ ^ 'TT^T) 

1.31(4) 

0.757122 

i 

c 

4.2(5)-10“'' 

5.2(4)-10“'' 


The partial decay width of the decay tt^ ^ 77 is thus 
given in our model by 

= Tra'^m-lo , (14) 


with 

T — 9(jJ7tp ' 9(jj^ 94>'Kp ' 94>'y V 9uj't^p ' 9(jj'^ V 9uj"'kp * 9(jj"^ 

(15) 

This ansatz leads to r^o = 6.6(3) • 10“^^ s for the tt^ ^ 
77 lifetime to be compared with the PDG El value r^o = 
8.3(6) • 10“^^ s. The uj' and uj" contributions are small but 
not negligible and the value of the tt^ lifetime obtained 
without them reads t^o{uj and 0 only) = 5.21(13) • 10“^^ s. 
Within two and a half standard deviations experiment 
and the model are in agreement {experiment — model = 
(1.7 ± 0.7) • 10“^^ s) and the remaining discrepancy can 
be attributed to the missing p-resonances or a more com¬ 
plicated Q^-dependence of the propagators. One has to 
remember that we make here extrapolation from the re¬ 
gion of cj, 0 , • • • resonances where the fit was performed 
to the region of the applicability of the chiral limit. Thus 
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we should recover here the chiral theory result for the 
TfO ^ lifetime (see [ 201 ) 


r{TT° 


77 ) 


9 

a m 0 

327r3/^ ’ 


(16) 


which gives r^o = 8.69 • 10“^^ s for f=132 MeV. Rephras¬ 
ing the same statement using the coupling constants one 
approximately expects 


1 

327t^P 




(17) 


where the dots correspond to terms neglected in modelling 
of the 37 r current. The discrepancy indicates that extend¬ 
ing the validity of our model beyond the description of 
the 37 r current one should probably take into account fur¬ 
ther contributions. However the 2.5 a discrepancy does 
not allow to draw any final conclusion. 

The amplitude of the process tt^ ^ 77 * for small val¬ 
ues of the four momentum of the off shell photon can be 
parameterized by a slope parameter a 

(1 T . (1^) 

Expanding the Breit-Wigner propagators one obtains 



QooTvp ' Quj'y 9(f>7Tp ’ 94)'y 


rpp Tml 

Quj'ttp ' 

Tm^, 


Tml 

9u}"7rp ’ 9uj"'y 

TmK 


(19) 


where T is defined in Eq. Numerically this gives a = 
1.74(2) GeV“^, or a • m^o = 0.0317(5), to be compared 
with the experimental value nn a • m^o = 0.032(4). 

The model also predicts the decay rates for the 
7 r^ 7 , uj 7 r ^7 and 0 ^ Tr^y: 


rpo_ 


ol c 


1 - 


m, 


2 13 


( 20 ) 




V ^77^7 


a 


V 


mlo 


1 3 


J 


2 2 
dVirp ’ 9p^ 5 


( 21 ) 


where V stands for uj or <p. The results are collected in 
Table 0] and compared with the respective experimental 
values. 

In the 7 r ^7 decay our prediction for this branch¬ 

ing ratio 5.2(2)T0“^ is in agreement with the 6.0(1.3)T0“^ 
of El within 1(7. As the only isospin violating effect in 
our model for this process is the charged-neutral pion 
mass difference, our prediction for the Br(/)^ ^ Tr^y) = 
5.2(2) • 10“^ is identical (within the errors) with the one 
for the neutral mode and is also in agreement with the 
data m] Br(p^ ^ ^r^y) = 4.5(5) • 10“^. 

For the (p 7 r ^7 decay the branching ratio Br(0 ^ 
7 r^ 7 ) = 0.99(12) • 10“^ is in agreement within 2cr with 
the value 1.23(12) • 10“^ from Ref. fH] . Our result for 


Br(cc; ^ 7 r^ 7 ) = 15.4(5)% overestimates however the mea¬ 
sured value ((8.921 q; 24)%) by factor 1.7. 

The extracted couplings determine also the cross sec¬ 
tion of the reaction e+e“ ^ 7 r^ 7 , 


r(e 




r%) = 




1 - 


mK 


9% 


3 V s 

I [^7u;7r/9 ■ Quj^ Qtp-Kp ■ 

~\~9u)'ttp ' 9oj''y T 9uj"ttp ■gu"^\ ■ BWp{s) 

‘ 9i07rp * * dcfi-Kp ‘ B^V 

'l~9uj'np 9 uj'7 ' -RHA;'(«s) T 9u)"7rp ' 9u)"'y ’ -RHA;"(«s)| 


( 22 ) 


with BWi = p^uj, (p, uoP uj" defined in Eq. 0. 

The comparison with existing data is shown 

in Fig. IHl where one standard deviation bands are given. 
Reasonable agreement is observed around the p resonance 
in contrast to the uj region. This is of course a reflection 
of the agreement and disagreement of the corresponding 
decay rates. A similar behaviour was also observed in Ref. 
m Their model predict a value for r{uj ^ Tr^y) well in 
accordance with the experimental value, but could not eas¬ 
ily reproduce r{uj ^ tt+tt^tt^). Further, both theoretical 
and experimental, studies of reactions e+e“ ^ tt+tt^tt^ 
and e+e“ ^ 7 r ^7 are thus required. 

The branching ratios of the cj, p and p decays to 
obtained within our model can be found in Table 0 They 
are in agreement with the PDG values El for a;, p within 
1-2 standard deviations. However the predicted value for 
Br(p^ ^ 7 r+ 7 r“ 7 r^) is more then two orders of magnitude 
smaller than the experimental value Br(p^ ^ tt+tt^tt^) 

1 • 10 “^, which is consistent with zero at two sigma level. 


Table 5. Branching ratios of the cj, p and p to tt^tt 7r° decays: 
model vs. experiment m 



model 

experiment 

Br(a; ^ tt^tt^tt^) 

95.1(22)% 

89.1(7)% 

Br(0 ^ TT^TT^TT^) 

14.5(22)% 

15.4(5)% 

Br(p ^ TT^TT^TT^) 

1.9(3)- 10"“ 

(1.01i^-||±0.34)- 10"^ 


4 Tests of the MC code 

In comparison to the previous versions of PHOKHARA 
the default random number generator was changed to the 
double precision version of RANLUX El written in G 
by Martin Liischer. A G-FORTRAN interface is provided 
with the PHOKHARA 5.0 distribution. 

To assure a technical precision of the code better than 
a fraction of one per mile a number of tests were per¬ 
formed for the new hadronic state in the generator. Among 
other tests, the initial state emission was tested against 
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e^e 7r°7 



V~s{GeV) 

e'^e~ 7r°7 



Data from Refs. |2ni (SND 2003) and |231 (SND 2000) are shown 
together with 1 a allowed bands for our model predictions. 


known analytical results, where all photon angles are in¬ 
tegrated, similarly to previously performed tests for other 
hadronic channels (see [7||Hlin0)- The independence of the 
result from the soft photon separation parameter was also 
checked. 


5 Summary and Conclusions 

The Monte Carlo event generator PHOKHARA has been 
extended to the three-pion mode. The model adopted for 


the hadronic form factor properly describes the currently 
available data for the cross section and the distributions. 
The ansatz is based on generalized vector dominance and 
is consistent with various other measurements like F (tt^ ^ 
77 ), the slope parameter of the tt^ ^ 77 * amplitude and 
radiative vector meson decays p 7 r^ 7 , (j) 7 r^ 7 , but is 

in conflict with uo Tr^y. 

The current version of the computer program (PHO¬ 
KHARA 5.0) is available at 

http://cern.ch/german.rodrigo/phokhara 
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